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ABSTRACT 



We map the full extent of a rich massive young cluster in the Cep 0B3b associ- 
^, . ation with the IRAC and MIPS instruments aboard the Spitzer Space Telescopqj 

tZt . and the ACIS instrument aboard the Chandra X-Ray Observatorja. At 700 pc, 

it is revealed to be the second nearest large ( > 1000 member), young (< 5 Myr) 
Ph, cluster known. In contrast to the nearest large cluster, the Orion Nebula Cluster, 

O ! Cep 0B3b is only lightly obscured and is mostly located in a large cavity carved 

'^ ! out of the surrounding molecular cloud. Our infrared and X-ray datasets, as 

^ ' well as visible photometry from the literature, are used to take a census of the 

young stars in Cep 0B3b. We find that the young stars within the cluster are 
>- ! concentrated in two sub-clusters; an eastern sub-cluster, near the Cep B molec- 

ular clump, and a western sub-cluster, near the Cep F molecular clump. Using 
O ! our census of young stars, we examine the fraction of young stars with infrared 

• ! excesses indicative of circumstellar disks. We create a map of the disk fraction 

m 
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throughout the cluster and find that it is spatially variable. Due to these spa- 
tial variations, the two sub-clusters exhibit substantially different average disk 
fractions from each other: 32% ± 4% and 50% ± 6%. We discuss whether the 
discrepant disk fractions are due to the photodestruction of disks by the high 
mass members of the cluster or whether they result from differences in the ages 
of the sub-clusters. We conclude that the discrepant disk fractions are most likely 
due to differences in the ages. 

Subject headings: pre-main sequence — stars: formation — infrared: stars 



1. INTRODUCTION 



Surveys of molecular clouds show that approximate ly half of low mass stars form in rich 



massive clusters w ith more than a thousand young stars (ICarpenteiil2000l : Ide Witt et al.ll2005 



Allen et al.ll2007| ). Most near by young clusters are of a small to moderate size, containing at 
most a few hundred members (JGutermuth et al.ll2009l ) . Consequently, much of what we know 
about star formation in rich, massive clusters comes fron a studying the neare st example: the 
Orion Nebula Cluster (ONC) at a distance of 414 pc (JMenten et al.l 120071 ). Here we use 
the term ONC to denote not just the dense, ~ 600 member cluster core found toward the 
brightest regions of the Orion Nebula, but rather the ex tended 8 pc x 4 p c region containing 
~ 2000 forming stars (Megeath et al. (in preparation); ICarpenterl ( 2000[ )). 



Ce pheus B3b was recognized as a distinct subgroup in the Cepheus 0B3 association by 



Blauuwl (Il964j ). Subsquent observations of the border between the Cepheus 0B3 molecular 
cloud and the Cep B3b subgroup detected a young group of > 400 low mass, pre-main 



sequence (PMS) stars (iMoreno-Corra 



Pozzo et al. 2003 



Getman et al 



et al.lll993l : Ijordi et al.lll996l : Mikami fc Ogural 12001 



20061 ) . Some of these s tars are embedd ed in the molecular 
clump Cep B (ISargentlll977l ) near the HII region Sh 155 (ISharplesslll959); how ever, most are 
lightly obscured and can be studied at visible wavelengths. IJordi et al.l (Il996[ ) established an 
age of 5.5 - 7.5 Myr by fitting theoretical isochrones to dereddened Fand B measurements 
for twenty- one stars in Cep OB 3b, consid erably older than the average age of the ONC 



< 2 Myr, iHillenbrandl 119971: iJeffriesI l2007h . Based on an analyis of the F vs. V-I color 



magnitude diagram, iMayne et al.l (l2007l ) proposed that Cep 0B3b has an age of ~3 Myr 
at an adopt ed distance of 8 50 pc . Using Chandra observations to identify PMS members of 
Cep 0B3b, ICetman et al.l (120091 ) determined extinctions to those stars with 2MASS near- 
IR color magnitude diagrams and fit evolutionary isochrones to a de-reddened V vs. V-I 
color-magnitude diagram of those stars. They deduced an age of 2-3 Myr for the part of the 



cluster near Cep B using an adopted distance of 725 pc. More recent work by iLittlefair et al. 
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( 120101 ) suggests an older age estimate (~4.5 Myr, compared to an age of 2 Myr for the ONC) 
with a distance of 580 ± 60 pc. Although the distance estimates in the literature vary from 
580-850 pc, we adopt a distance of 700 pc; this is the maser parallax-deterni ined distance for 



Cepheus A which lies in the same molecular cloud (JMoscadelli et al.ll2008[ ). In all of these 
studies, Cep 0B3b appears to be older than the ONC. 

As part of a survey of the Cepheus 0B3 molecular cloud we have mapped the Cep 0B3b 
cluster with the IRAC and MIPS instruments aboard Spitzer. When the Spitzer observations 
revealed the full extent of the cluster, we then mapped the cluster with the ACIS instrument 
aboard Chandra. This is the first mid-infrared and X-ray census of young, low mass stars 
in the entire Cep 0B3b cluster. Using these data, we map the distribution of young stars 
with enhanced infrared (due to dusty circumstellar material) or X-ray (due to increased 
chromospheric activity) emission. These data show that Cep 0B3b is one of the largest 
known clusters within 1 kpc of the Sun. We propose that the cluster may be a rare, nearby 
analog of the ONC at a more advanced stage of evolution. In this paper, we present a census 
of young stars, both with and without disks, as well as protostars in the Cep 0B3b cluster. 
Using this census, we show that the fraction of young stars with disks varies within the 
cluster. Finally, we discuss reasons for this variation. 



2. OBSERVATIONS & DATA ANALYSIS 
2.1. Mid Infrared Data 

Mid-infr ared (MIR) observations were obtained with the Spitzer Infrared Array Cam- 



era (IRAC) (IFazio et all \2004] and the Multi-band Imaging Photometer System (MIPS) 
(JRieke et al.ll2004l ) as part of a Spitzer Cycle-2 GO survey (PID 20403) of the entire Cepheus 
0B3 molecular cloud. The molecular cloud was mapped by IRAC at two different epochs; 
the first on 2005 December 25 and the second 2006 August 8. The map was further ex- 
tended in two epochs (2007 August 17 and 2008 January 31) to cover the entire Cep 0B3b 
cluster as part of the Spitzer program PID 40147. Each field was observed four times in 12 s 
High Dynamic Range (HDR) mode. A total integration time of 41.6 s was obtained at each 
wavelength for each field. 

The Spitzer /MIPS 24 /im images were obtained on 2007 February 25-26 and 2007 Octo- 
ber 27 as part of the programs described above (PID 20403 and 40147), using scan maps at a 
medium scan rate with an offset of 302" between scan legs. Complete coverage was obtained 
only for the 24 /xm band; consequently, we only present results from this MIPS band. This 
paper focuses on a 49' x 33' region on the north-west edge of the large Cep 0B3 molecular 
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cloud which encompasses Cep 0B3b. 

Mosaics of the IRAC data were constructed from the basic cahbrated data (BCD) prod- 
ucts provided by the Spitzer Science Center using their standard data processing pipehne 
versions S18.5 and 18.7. Post-BCD reduction of the Spitzer/IRAC da ta was executed using a 
software package called "Cluster Grinder" (see lGutermuth et al.ll2009l ). developed specifically 
to work on young clusters. The package includes modules for bright source artifact correc- 
tion, cosmic ray remova l, mosaicking, and the visualization and photometry tool "Photvis" 
( JGutermuth et al.ll2008l ) for point source extraction and aperture photometry. Radii of 2". 4, 
2". 4, and 7". 2 were used for the aperture, inner, and outer sky ann ular li mits, respectively 



(choic es of aperture and annuli are justified in lAUen et al.l (|2008[ ) and iGutermuth et al. 



The instrumental magnitudes were derived directly from the BCD frames converted 
into units of DN. For the 10.4 second frames, the adopted zero-magnitudes were 22.00, 



21.24 , 19.04 and 19.43 mag for the 3.6, 4.5, 5.8 and 8 fim. bands, respectively ( JReach et al. 



20051 ). These included aperture corrections from a 2 pixel to 10 pixel aperture of 1.21, 1.23, 
1.58, and 1.58, respectively, as tabulated in the IRAC Data Handbook. The Spitzer band 
magnitude labels are denoted by [wavelength in microns]: for example, [3.6] = 15.16 means 
that the 3.6 /xm magnitude is 15.16. Ninety percent differential completeness limits of the 
final mosaics are approxima tely [3.6] = 15.16, [4.5] = 15.18, [5.8] = 13.71, and [8.0] = 11.72 



see 



Gutermuth et al.l ( 120051 ) for the estimation technique). 



PID 20403 MIPS da ta were processed by the MIPS instrument teams' Data Analysis 
Tool ( JGordon et al.ll2005l ). which calibrates the data, corrects distortions, and rejects cosmic 
rays during the co-adding and mosaicking of individual frames. The MIPS 24 fira mosaic 
created from the four scan maps has a pixel size of 1".245 and units of DN s~^ pix~^. The 
MIPS field was extended in the PID 40147. For this second observation, MIPS BCD data 
from the S16.1 pipeline were mosaicked using the "Cluster Grinder" package. A pixel size 
of 2.547" was used for this mosaic. Point sources were identified in the Spitzer/MIPS data 
using Photvis; aperture photometry was then obtained using 7". 4, 14". 7, and 24". 5 for the 
aperture, inner, and outer sky annular limits. These magnitudes were then used as the initial 
guesses for point spread function (PSF) fitting photomet ry of the identified sources using 
the IDL implementation of DAOPHOT ( jLandsmaru Il995l ) : this vers ion of DAOPHOT wa s 
modified to mask out saturated pixels in partially saturated PSFs (JKryukova et al.ll2012l ). 
Two point spread functions, one to fit unsaturated sources and a larger one to be used to fit 
partially saturated sources, were created from fifteen sample stars located across the MIPS 
field. The PSF generated for the unsaturated sources had a size of 62" x 62", while the 
PSF generated for the saturated sources was 87" x 87" in size. The fitting radii were 2.5" 
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and 12.45" for the unsaturated and saturated PSFs, respectively. The large fitting radius 
for the saturated PSF was needed to fit the wings of the PSF after the central, saturated 
pixels were excluded from the fit. The MIPS photometry were calibra ted using a zero-point 



flux density of 7.17 Jy, referenced to Vega (JEngelbracht et al.l 120071 ). The ninety percent 



differential completeness limit is [24] = 7.68. A three-color rendering of the Spitzer data can 
be seen in the top pane of Fig. [H 



2.2. X-ray Data 

X-ray data were obtained with the ACIS instrument aboard the Chandra X-ray Obser- 
vatory. Two 17' X 17' fields were imaged covering both sub-clusters. Three images of the 
field covering part of the eastern sub-cluster were obtained with an aim point of 22:55:47.5, 
+62:38:10.3 (PIDs 9920, 10809 and 10810, with roll angles of 48°, 35° and 72° respectively), 
and three images of the field covering the western sub-cluster were obtained with an aim 
point of 22:53:31.7, +62:35:33.4 (PIDs 9919, 10811 and 10812, with roll angles of 68°, 63° 
and 63 respectively). Each image has an integration tim e of 28 ks. The rema inder of the 



eastern subcluster was previously imaged with Chandra by lGetman et al.l (120061 ) (PID 3502, 
aim point of 22:56:49.4, +62:39:55.6 and roll angle of 7.9 ). The fields are outlined in the 
bottom pane of Fig. [H 

These data have all been reduced in a similar manner using the Chandra Interactive 
Analysis of Observations (CIAO) software tools. The data were processed through the 
standard CIAO pipeline at the Chandra X-ray Center, using their software version DS7.6. 
This version of the pipeline automatically incorporated a noise correction for low energy 
events. This filter can remove good events from the cores of bright point sources, resulting 
in an underestimation of the X-ray flux. In this case, the count rate did not exceed 0.02 
counts per second for any source thus, event loss was not a concern. Background was nominal 
and non- variable. 

As the goal of the observations was the correlation of X-ray and infrared sources, point 
source detection was aggressive. First the CIAO script fluximage was used. The fiuximage 
script creates exposure-corrected images of ACIS observations as well as the corresponding 
aspect and instrument maps which prepare the data for point source detection. The script 
is especially efficient in cases in which multiple observations are used. Sources were detected 
using the Wavdetect routine. Wavdetect correlates the image with point spread functions 
(wavelets) of different scales and then searches the results for significant correlations. In 
this case, scales of 1 to 16 pixels were search on with logarithmic spacing. Significances as 
low as 1 sigma were returned. To ensure detection, we require a corresponding detection 
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in the infrared to confirm tlie X-ray nature of the source. As discussed in § 13. 2[ the vast 
majority of X-ray sources with visible and infrared counterparts are found toward the Cep 
0B3b isochrone. THis shows that there are few spurious X-ray sources in our sample. 



2.3. Near Infrared and Visible Data 

Sources detected by Spitzer and Chandra were matched with c ounterparts from th e 



2MASS Point Source Catalog (PSC) and the BVI band catalog of iMavne et all (120071 ). 
A matching radius of 2" was used in all cases. The addition of these catalogs extends 
the wavelength coverage for these objects into the near-infrared (~1 fim) and the visible 
(~0.5/im). A catalog combining the Spitzer, 2MASS and Chandra data will be called the 
X-ray sample and a catalog combing the Spitzer, 2MASS and VBI data will be called the 
visible sample. 



3. THE IDENTIFICATION OF YOUNG STELLAR OBJECTS 
3.1. IR Identification of Young Stars with Circumstellar Material 



Gutermuth et al.l (120091 ) presented a new empirical three-phase scheme to identify young 
stellar objects (YSOs) with infrared excess (i.e. protostars with dusty infalling envelopes and 
young pre-main sequence stars with dusty disks) which we adopt here. Those classifications 
are illustrated in the IRAC/MIPS color-color and color-magnitude diagrams shown in Fig. [2l 

In phase 1, all four IRAC bands with photometric uncertainties a < 0.2 mag are used 
to separate sources with infrared excesses due to disks and infalling envelopes (Class II and 
Class O/I, respectively) from pure photospheres (Class Ill/field stars) based on [3.6] — [5.8] 
and [4.5] — [8] colors. Various non-stellar contaminating sources (e.g. PAH and AGN extra- 
galactic sources) also s how infrared exce s ses; t hey must be eliminated by several criteria 



described in detail by JGutermuth et al.l (12008! ) . Appendix A. Galaxies with strong PAH 



emission are eliminated by their colors. AGNs mimic the colors of young stellar objects and 
must be elimiminated by their faintness, as determined by their location in the [4.5] vs. [4.5] 
- [8] color- magnitude diagram (these cuts are shown in Fig. [2]). Unresolved shocked emission 
blobs arising in YSO outflows are identified by their strong 4.5 /xm emission. All sources 
classified as contaminants were re-examined in the third phase of our identification process for 
excess emission at 24 fun since some protostars might be masquerading as contaminants. The 
remaining infrared excess sources which are not removed as contaminants can be classified by 
[4.5] - [5.8] color as Class O/I objects (protostars) or Class II objects (pre-main sequence star 
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with disk). The sources that do not show IR-excesses are consistent with pure photospheric 
emssion between 1 and 8 fim. These include Class III (diskless pre-main sequence star), field 
stars, as well as transition disk objects (stars with disks exhibiting AU-scale inner holes) 
which will be identified in later phases. The distinct mid-IR colors of Class O/I objects. 
Class II objects and pure photospheres are clearly evident in Fig. |2l 

Phase 2 is applied to those sources which lack IRAC detections at 5.8 or 8.0 /im, but 
which have high quality 2MASS detections (cr < 0.1) in at least H and Kg bands. Line 
of sight extinction to each of these sources is estimated using 2MASS colors (if J, H, Kg 
detections are all av ailable) by de - redde ning the sources that lie above the classical T-Tauri 
star (CTTS) locus of iMever et al.l (Il997l ) or an adopted standard dwarf-star color {{J—H)q = 
0.6, iBessel &: BrettI Il988l ) in J — if and H — Kg color space. The sources that lie on or 
below these fiducial values are assumed to have no measurable reddening. If high-quahty J 
detections are unavailable, the line of sight extinction is estir nated by deredden ing the source 
to the CTTS loc us in H — Kg an d [3.6] — [4.5] color space (see lAUen et al.ll2008l ) and using the 
reddening law of Flaherty et al.l (120071 ): this law was derived for star forming regions similar 
to Cep 0B3b. Sources are then classified as either Class Ill/field objects. Class II objects 
or Class 0/1 objects based on their position in dereddened {K — [3.6])o and ([3.6] — [4.5])o 
color space ( lAUen et al.ll2008l ). It must be noted that this phase applies primarily to lower- 
mass objects with disks, which are not detected at the longer IRAC wavelengths due to 
the increased background levels at those wavelengths. Due to the rising spectral energy 
distribution (SED) of the protostars with envelopes, we expect most of them to be detected 
at the longer IRAC wavelengths and thus, we expect this phase to only add a few, if any, 
protostars to our catalog. As expected, all the protostars in our catalog are classified during 
phase 1. 

In addition, there are a handful of objects with Near-IR-|-IRAC colors of Class Ills, but 
which exhibit a 24 fxm excess (e.g. ([5.8] - [24]) > 1.5 and ([8] - [24]) > 1.5). These objects 
are classified as transition dis k objects, pre-main seq uence stars surrounded by disks that 
contain AU-scale inner holes ( JMuzeroUe et al. Il2010l ). Due to the inner hole, these disks 
typically emit weakly at wavelengths < 8 /im. 

Phase 3 is applied to all sources classified in phase 1 with MIPS 24 fim detections with 
0" < 0.2 mag. This phase re-examines all sources including contaminants and identifies Class 
0/1 objects and transition disk objects on the basis of their 24 /im photometry. In this phase. 
Class O/I objects are defined as sources with [5.8] - [24]) > 4. The [3.6] - [4.5] vs. [4.5] - [24] 
color-color diagram (Fig. [2]) clearly shows the distinct colors of our identified Class O/I and 
Class II sources. 



In total, there are 1135 objects in the Near-IR -|- IRAC dataset with evidence for dusty 



envelopes or disks: these objects will be used in the following analysis. Of these, 56 are 
classified as protostars (Class O/I), 1014 as pre-main sequence stars with disks (Class II), 
and 65 are transition disk objects with AU-scale disk holes. Tables [H [21 [Hand ?? list the 
sources identified as Class O/I, Class II, transition disk and Class III, respectively. A full 
source list is included with the electronic version of this paper. 



3.2. X-ray Identification of Young Stars with Enhanced Coronal Activity 

An estimate of the total number of young stars both with and without disks requires 
a method for identifying and counting diskless Class III objects. These pure photospheric 
sources cannot be distinguished from background field stars by their near to mid-IR colors. 
In the following sub-sections, we use two separate methods to quantify the number of diskless 
Class III objects. The first method directly identifies potential Class III objects from their 
X-ray flux density, and the second method distinguishes potential Class III objects using 
their visible magnitudes and colors. 

Young stars have enhanced X-ray luminosities by orders of magnitude wit h respect to the 
main sequence (MS) in large part due to increased chromospheric activity (JFeigelson et al. 



20071 ). By selecting objects that have both photospheric colors in the infrared as well as X- 
ray emission, we can identify diskless young members of the cluster. To make this selection, 
we only need the spatial positions of the X-ray sources. The analysis of the flux densities 
and other X-ray properties of these sources will be discussed in a future paper. 

We identify as X-ray members of Cep 0B3b all sources that have a Chandra X-ray 
detection, and have been classified as a Class 0/I/II/photospheric source using the Spitzer 
and 2MASS three-phase classification scheme. However, it is possible that coronal activity 
from another type of X-ray active field star, such as a flaring dMe star (the X-ray luminosity 



of dMe stars can reach up to log L^ ~ 30 erg s ^ (JMarino et al.lll998l )) could lead to a 



misidentiflcation of a fleld star as a young star. In comparis on, the typical young star has 



an X-ray luminosity in the range of 28 < log L^. < 32 erg s~^ (JFeigelson et al.l 120071 )). Thus, 
the X-ray active fleld stars may contaminate the Chandra sample and be misclassifled as a 
Class III source. To estimate the potential magnitude of the background contamination, we 
examine the location of the X-ray sources in V and V — I color-magnitude space (right pane 
Fig.E]). We consider the sub-set of sources that have detections in all three data sets Spitzer, 
Chandra, and BVI. This sample contains: Class O/I, 136 Class II, 9 transition disk and 
329 Class III sources. We assume that the potential Class III sources that lie below the zero 
age main sequence (ZAMS) in Fig. |3] are all background contaminants (ther e are only two 



Class II sources in that region, they are the sources with V — I < 2). Further iGetman et al. 
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( 120061 ) estimate that 13 sources in their 17' x 17' of Cep B are foreground stellar sources. If 
we apply this number to out full Chandra field of view, which is 729 arcminute^, we estimate 
that there will be 33 additional contaminants. Combining these two estimates yield 11% 
(53/499) of the sources as contaminants. Considering the full Chandra field of view of 729 
arcminute^, we have an X-ray contamination rate of 0.07 stars arcminute^. 

The X-ray sample includes: 5 Class O/I, 233 Class II, 16 transition disk, and 499 Class 
III sources, of which 53 may be contaminating X-ray active field stars. For reference, within 
the region imaged by Chandra there are 36 Class O/I, 810 Class II, 29 transition disk and 
3213 photospheric sources (Class III or field stars) as classified via Spitzer and 2MASS colors 
alone. Thus, we are detecting 29%(233/801) of the Class II sources and 55%(16/29) of the 
transition disk sources with Chandra. 



3.3. Visible Wavelength Selection of Young Stars 

We can also search for diskless young stars by utilizing the intrinsic magnitudes and 
colors of young stars to separate them from field stars in l^ vs. V — I space (Fig. [3]). 
Young stars are still contracting to the MS, thus for a given temperature, young stars have 
a greater luminosity than field stars at an equal or greater distance. The young stars occupy 
an isochrone (or bands of isochrones) that lies above the MS, and the extinction vector is 



approximately parallel to this isochrone. Using the V and / band data of lMayne et al.l (120071 ) 
for stars with V and V — I photometric errors less than 0.1 mag, we empirically determine 
the YSO locus in V and V — I space. To determine the isochrone of young stars in the 
cluster, the locations of members identified by Spitzer and Chandra are plotted in V and 

V — I space (Fig. [2]). The young stars are then divided into 0.75 magnitude bins between V 
magnitudes of 15 and 24, and a two step iterative process is used to determine the median 

V — I color of the known YSO population in each bin. First, the median color of the bin is 
calculated using all sources in the bin. Next outliers greater than 3a from the median are 
rejected and a new median and standard deviation are calculated. The width of the YSO 
isochrone in color space is then taken to be between — 1.25cr and 2cr around the median color 
of the bin. The median color of the YSO isochrone is fit with a third order polynomial, 
as are the outer bounds of the YSO isochrone; all stars between the two outer bounds are 
considered to be likely members. The lower bound of the isochrone is cropped closer to the 
median value because the density of sources quickly increases in the direction perpendicular 
to the isochrone in the direction toward the ZAMS. In particular, foreground dwarfs and 
background giants may masquerade as YSOs by appearing more luminous for a given color. 
The outer bounds are shown in Fig [31 Hereafter, this selection of sources will be denoted by 
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the visible sample. 

We can estimate an upper bound to the possible non-member stellar contamination 
(bot h foreground and b ackground) by considering the V and / photometry from field 3a 



from IPozzo et al.l (120031 ). hereafter called the Pozzo field. This is a 12' x 12' field centered 
on RA: 22:50:21.14 and Dec: 62:46:44.88, on the north-west of the molecular clump Cep 
F. We choose this region due to its location off the cluster and because it has the least 
number of Spitzer identified excess sources; only four. If we assume that these four sources 
are the only sources in this field that are members of Cep 0B3b, then we can apply the 
same V and V — I selection criteria and consider all remaining sources that occupy the YSO 
dominated region of color-magnitude space as contaminants. There are 98 sources that fit 
this criteria (102 sources in the YSO dominated region, minus the 4 known excess sources), 
for a contamination rate of 0.6 stars arcminute"^. 

Another way to estimate the number of non-member stellar contaminants in the V and 
V — I space is to extrapolate from the number of sources in the region of V and V — I 
space below the YSO isochrone into the YSO isochrone. Assuming that there is a constant 
stellar density along t he line of sight, then the n umber of stars of a given magnitude, m, is 



approximately 10"°'" (iBinney &: Merrifeldlll998l ). We approximate the YSO isochrone with 



a linear strip in V and V — I space that is two magnitudes wide (bounded by the lines 

V = 2.9(\^ ~ ^) + 9.4 and V = 2.9(y — I) + 11.4), and consider the sources found in the one 
magnitude linear strip below the YSO isochrone,bounded by the lines V = 2.9{V — I) + 11.4, 

V = 2.9(1^ ~ -^) + 12.4, V = 14.5 and V = 22.5 to be non-member field stars. Then the ratio 
of non-member field stars along the line of sight within the YSO isochrone to the number of 
stars within the field star isochrone is just 

AT {'^ ^C\^■^"^^\m 

-''contaminants Jq >-ii/(/ r\ Q-\ A /i \ 

A^fieid ~ J^ 100-6'"dm ~ ■ 

Since there are 1211 sources in the one magnitude field star strip, we estimate 381 contami- 
nants in the YSO isochrone. Dividing by the area covered by the visible data (4 x 11' x 22') 
leads to a contamination rate of 0.4 stars arcminute"^. This figure is similar to, but slightly 
lower than, the density of contamination derived from the Pozzo field. Because the Pozzo 
field may contain some young sources due to its proximity to the molecular cloud, we adopt 
0.4 stars arcminute"^ as the visible selection contamination rate. 

The visible sample includes: 402 Class II, 23 transition disk, and 1270 Class III sources, 
of which 381 may be contaminants. If we restrict the sample to only the sources that lie 
within the region with a field of view in common with the Chandra images (we will call this 
region the overlap region), we find 372 Class II, 17 transition disk, and 914 Class III sources, 
with an estimated 232 contaminants. 
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4. THE SIZE AND STRUCTURE OF CEP OB3B 

The entire population of identified young stars (see Fig. H^) is found in a ^-^ 7pc x lOpc 
region. The Class I and II objects are concentrated in two sub-clusters (hereafter, the eastern 
and western sub-clusters), surrounded by an extended halo of these objects. Furthermore, 
the two sub-clusters are located primarily in a cavity with surfaces defined by bright PAH 
emission; this cavity was probably evacuated by the OB stars in the cluster (see Fig. [1^). 
The eastern sub-cluster has a Class II surface density peaking at about 260 pc~^ and the 
western sub-cluster has a peak Class II surface density of about 310 pc^^ (see Fig. |Dd). The 
sub-clusters are separated and surrounded by a diffuse 'halo' of young stars. We define the 
halo as the region with a Class II surface density above 10 pc~^. The total area of the cluster, 
including the diffuse halo is about 70 pc^. There are 837 Class II sources within this area 
with an average surface density of about 12 pc~^. 

If we include the Class III YSOs determined via Spitzer and X-ray detections, the peak 
YSO densities reach ~ 340 pc~^ and ~ 680 pc~^ for the eastern and western sub-clusters, 
respectively. This sample has 1317 sources within the 70 pc^ region defined above as the 
total area of the cluster, for an average surface density, including the halo, of about 19 pc~^. 
If we consider the X-ray contamination rate of 0.07 stars arcminute^ these surface densities 
will be reduced by about 10 pc^^. 

If we consider the Spitzer plus visible sample, the peak surface densities are ~ 540 
pc~^ and ~ 730 pc~^ for the eastern and western sub-clusters, respectively, and the average 
surface density including the halo is ~ 30 pc^^ (2112 sources / 70 pc^). Including the visible 
sample contamination rate of 0.4 stars arcminute"^, these values will be reduced by about 
50 pc-2. 

An extinction map (see Fig. HJ top pane) was derived using the H — Kg color s of back- 



groun d stars from the 2MASS (Two-Micron All Sky Surv ey) Point Source Catalog (ISkrutskie 



20061 ). The method utilized is outlined in detail in §3.1 of ICutermuth et a. 



porat es elements of the NICE (JLada et al.lll994j ) and NICER algorithms (JLombardi fc Alves 



(2005) and incor- 



20011 ). Briefly, the line-of-sight extinction to each point on a 31" grid is estimated by evaluat- 
ing the mean and standard deviation of the H — Kg colors of the 20 nearest stars to the grid 
point, and an iterative outlier H — Kg color algorithm was used to reject colors more than 
3cr from the mean until the mean converged. The extinction map was generated assuming 
an average H — Kg color of 0.2 for background stars, A^^ = I.S2{Eh-ks), and Ay = 9.9Ak^ 



(JRieke fc Lebofskylll985l ). Typical resolution is 2', and the range of Ay is - 12.6 mag in 



this map. 

The sub-clusters are associated with two distinct molecular clumps: Cep B to the east 
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and Cep F to the south-west using the nomenclature of ISargentI (Il977l ). As noted, a large 
percentage of the pre-MS population is in the low extinction cavity: 89%(965/1079) of Class 
II and transition disk objects are projected on a region with an Ay < 5 based on the derived 
extinction map. This observed configuration is different from that found in the ONC While 
the ONC is in front of and partially embedded in the molecular cloud containing BN/KL, 
Cep 0B3b lies mostly adjacent to two active molecular clumps. This is advantageous because 
the detectability of low-mass pre-main sequence stars in Cep 0B3b by Spitzer will be less 
affected by extinction and nebulosity than in the ONC. 

The 56 protostars (i.e. sources classified as Class O/I) are primarily found in the vicin- 
ity of the Cep B and Cep F molecular clumps. About 41%(23/56) are associated with high 
extinction, Ay > 5, based on the derived extinction map. Tis shows that star formation is 
continuing in these clumps. The remaining 33 protostars apparently are not associated with 
high extinction, although this may be a short-coming in our extinction map because of its low 
spatial resolution, and the consequent bias against detecting compact hig h extinction struc 



tures. On the other hand, some of these protostars may be edge-on disks; iGutermuth et al. 



(120091 ) estimated that the number of edge-on disks is at most 3.5% of the total number of 
Class II objects. In Cep 0B3b this upper limit corresponds to ~35 edge-on disk systems. 
Since this number is comparable to the number of protostars in Cep 0B3b that are not 
associated with high extinction, most of these YSOs classified as protostars may be edge-on 
disk systems. 



5. THE DISK FRACTION IN CEP OB3b AND ITS SUB-CLUSTERS 

5.1. Disk Fraction in the Sub-Clusters 

To estimate the disk fraction we require a census of the cluster members with and 
without disks. While the number of members with disks is reliably determined from our 
sample of IR excess sources, to determine the number of diskless members we require a 
method to distinguish between diskless Class III and field stars seen in projection toward 
the cluster. This has often been done by using measurements of the source density in off- 



cluster regions to estimate the number of field stars (see for instance ICarpenterl (J2000l ) and 



Allen et al.l ( 120071 )). However, these methods require assumptions about the homogeneity 
of the background distribution of stars and line of sight extinction that may not hold in a 
large, complex OB association with multiple epochs of star formation. Thus, we use the 
direct identification of Class III members using either their X-ray flux or their position on 
the V^ vs. V — I diagram, as described in section 13.31 
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Since disk sources are selected because of their increased flux in the infrared, sources 
with disks can be detected to fainter magnitudes, and higher extinctions, than purely pho- 
tospheric sources. To reduce the impact of this bias, we apply a near-IR magnitude and 
color requirement based on the position of the source on the J vs. J — H color-magnitude 
diagram (Fig. [5]). The sources that remain are those that are brighter than a de- reddened 
J-band magnitude of 14.75. This magnitude an d reddening vector c orresponds to a 0.2 M© 



star at a distance of 700 pc using the models of iBaraffe et al.l (119981 ). Further, these sources 
must have colors that indicate an extinction of less than 6 Ay- These limits were set based 
on a visual inspection of the positions of the sources in Fig. [51 

After applying these magnitude and extinction cuts, we are left with 51 Class O/I, 727 
Class II and 56 transition disk sources from the Spitzer and 2MASS sample. For the X-ray 
sample 2 Class O/I, 207 Class II, 16 transition disk and 451 Class III sources remain, and 
for the visible sample 345 Class II, 22 transition disk, and 1139 Class III sources remain. 

We calculate the disk fraction for the entire cluster and for two 5' radius regions cen- 
tered on the eastern (RA: 22:56:49.44, Dec: 62:39:55.54) and western (RA: 22:53:50.83, Dec: 
62:35:36.37) sub-clusters, shown as red circles in Fig. S^. The disk fraction for the entire 
X-ray sample is 33% ± 2% (207 Class II, 16 transition disk and 451 Class III). The eastern 
sub-cluster has a disk fraction of 32% ± 4% (62 Class II, 2 transition disk and 137 Class 
III), while for the western sub-cluster, the disk fraction is 50% ± 6% (64 Class II, 7 transi- 
tion disk and 71 Class III). The disk fraction for the visible sample is 24% ±1% (345 Class 
II, 22 transition disk and 1139 Class III). This is significantly lower that the disk fraction 
estimated using the X-ray sample. The visible data used in this analysis, however, extend 
to the north of the region imaged by Chandra. If we restrict our estimate to only those 
sources in the visible sample located in the region imaged by Chandra, the disk fraction is 
29% ±2% (321 Class II, 17 transition disk and 811 Class III). This is consistent with the disk 
fraction estimated with the X-ray sample. Further, we find disk fractions of 34% ±3% - (i.e., 
103 Class II, transition disk and 199 Class III) in the eastern sub-cluster and 44% ± 5% 
- (i.e., 71 Class II, 4 transition disk and 96 Class III) in the western sub-cluster. Note that 
the fraction of sources with disks is similar regardless of the method used to identify the 
diskless Class III sources. These results are summarized in Table where we tabulate the 
disk fractions with and without the J vs. J — H cuts. 

Although we have minimized the contamination of the Class III sample by requiring 
either an X-ray detection or V vs V-I colors coincident with the CepOB3b isochrone, the 
residual contamination is the primary systematic uncertainty in our disk fraction measure- 
ment. In Table O we correct for the background contamination using the 0.4 arcminute"^ 
contamination rate from § 3. Since we do not know the J and if -band magnitudes of the 
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contamination, we do not apply the J vs J — H cut and background subtraction simulta- 
neously. The correction for contamination raises the disk fraction because the background 
subtraction is only applied to the Class III objects. For example, the background subtracted 
disk fractions for the overlap region increases to 36% ± 2% for both the X-ray and visible 
sample. Although the sub-clusters subtend a smaller region of the sky and hence have smaller 
rates of contamination; the subtraction of contamination can accentuate the difference in the 
disk fractions. In particular, if we consider the visible selection contamination rate of 0.4 
stars arcminute^, this is equivalent to 31 contaminants in each sub-cluster. After subtracting 
these from the Class III count, the disk fractions increase to 39% ± 4% - (i.e., 121 Class II, 
transition disk and 191 Class III) in the eastern sub-cluster and 53% ± 6%53 - ( i.e., 83 
Class II, 4 transition disk and 75 Class III) in the western sub-cluster. If we subtract out 
the potential contamination in the X-ray data (0.07 arcminute"^) then this corresponds to a 
contamination rate of 6 contaminating sources in each sub cluster, resulting in disk fractions 
of 34% ± 4% in the eastern sub-cluster (72 class II, 2 transition disk and 145 class III) and 
52% ± 6% in the western sub-cluster (73 class II, 7 transition disk and 73 class III), see 
Table El 

We place the disk fraction estimates for Cep 0B3b and its two sub-clusters within the 
context of other star fo rming regions by placing them on a plot similar to Figure 14 from 



Hernandez et al.l ( 120081 ). as seen in Fig. |6l We adopt the disk fractions for the X-ray sample 
with the J vs. J — H cut (32% ± 4% in the eastern sub-cluster and 50% ± 6% in the western 
sub-cluster); the systematic uncertainties due to Class III selection and contamination are 
similar to the formal uncertainties displayed in the figure. The age is assumed to be between 
3 and 5 Myr, represented by the horizontal uncertainty bars on the plot. 



5.2. Disk Fraction as a Function of Dereddened V-I Color 

Spitzer observations of young stellar clusters suggest that the disk fraction may be 
mass dependent, with the gas rich disks around more the massive stars disappearing first 



( jLada et al.ll2006t [Hernandez et al.l 120071 12009| ). To search for a similar trend in our data 



we use the dereddened V — I color as a proxy for the mass, where sources are dereddened 
with the line of sight extinction estimated using the method described in § 13.11 . To directly 
compare the visible and X-ray samples, we consider those members detected in the V and 
/-bands in the field of view common to both the visible and X-ray samples. We bin the 
sources into three bins of blue, intermediate and red colors, with roughly 380 sources in 
each bin. The dereddened V — I bounds for these three bins are: -0.5 to 1.05 mag, 1.05 
to 1.72 mag and 1.72 to 4.35 mag, respectively. We use the blue, intermediate and red 
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bins as proxies for high, intermediate and low-mass young stars. To roughly estimate the 
correspond ing mass ranges for to these color bins, we compare these dereddened colors to 
Table 5 in iKenyon &: Hartmann I (Il995l ) . These colors correspond to spectral t ype ranges 



of gre ater than K3, K3 to MO and MO to M6, respectively. Using the tracks of ISiess et al. 



( I2OOOI ) and an age of 3 Myr, these correspond to masses of greater than 1.7 Mq, 1.7 to 0.6 M© 
and 0.6 to 0.2 Mq. Table [6] shows the results for the field of view common to the X-ray and 
visible light observations and for the sub-clusters, respectively. 

We would also like to estimate the potential number of non-member stellar contaminants 
in each of these color bins. To do so we consider the same field star reference strip as defined in 
§ 13.31 However, to determine the color bins above, we depend on the line of sight extinction 
estimates to each of the sources. We do not know what these will be for the potential 
contaminants. Instead we re-redden the bounds for the color bins using the median Ay of 
all the sources considered to create the color bins. The median Ay of these sources is 3.3 
mag, and the resulting V ~ I color bounds are 1.22 to 2.77 mag, 2.77 to 3.44 mag and 3.44 
to 6.07 mag for the blue, intermediate and red bins, respectively. Since there are 739, 340 
and 115 sources in the blue, intermediate and red sections of the field star strip, we estimate 
232, 107 and 36 contaminants using Equation 1, and a contamination rate of 0.24, 0.11 and 
0.04 stars arcminute"^, respectively. 

Table [6] shows the disk fraction as a function of dereddened V — I color bin. For the 
common field of view, the disk fraction is the same for all three bins of the X-ray sample, 
roughly 30% ±5%, whereas for the visible sample, the intermediate bin (31% ±3% for 1.7 to 
0.6 Mq) has a slightly lower disk fraction than the reddest bin (38% ±3% for 0.6 to 0.2 Mq), 
and the blue bin (19% ± 2% for > 1.7 Mq) shows a substantially lower disk fraction than 
the intermediate bin. Note that if we subtract the potential contaminants from the visible 
sample, the large discrepancy between the visible sample disk fraction (30% ± 4%) and the 
X-ray sample disk (30% ± 5%) in the blue bin disappears, while for the red bin the disk 
fractions differ by about 1.7 a (41% ± 3% for the visible sample and 31% ± 5% for the X-ray 
sample). Examination of the right pane of Figure |3] shows that the discrepancy in the red 
bin could be due to a lack of X-ray sensitivity to the reddest, lowest mass members. 

Looking at the sub-clusters in Table O we find that the disk fraction does not vary 
greatly between the three bins for the western sub-cluster. In contrast, the disk fraction 
in the blue bin (19% ± 4% > 1.7 M© ) for the eastern sub-cluster is half that of the disk 
fraction in the intermediate bin (40% ± 6% 1.7 to 0.6 Mq). If we examine the difference 
in the disk fraction between the two sub-clusters, this difference is greater that 2cr for both 
the visible sample and the X-ray sample in the blue bin (> 1.7 Mq). The disk fraction can 
be found in Table [61 In contrast, for the intermediate bin (1.7-0.6 Mq) the disk fractions 
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of the two sub-clusters are within la for the visible sample and 1.2 a for the X-ray sample. 
For reddest bin (0.6 to 0.2 Mq), the difference in disk fractions is about 1.3a for the X-ray 
sample and 1.2cr for the visible sample. Subtracting the potential contaminants from the 
visible sample increases the difference in the disk fraction between the two sub-clusters in the 
bluest bin. The intermediate bin still has no difference, and the red bin has a difference of 
l.lcr. Thus, the lowest disk fractions and the largest difference in the disk fraction between 
the sub-clusters is in the bluest V — I bin with the highest mass stars. As we will discuss in 
§ 16.31 the observed difference for the bluest stars may be due to a systematically older age 
in the eastern sub-cluster. 

Broad-band magnitude, and color, determinations of stellar mass are highly uncertain 
due to the number of non-mass related factors that could affect the magnitude, and color, of 
a young source, such as age, accretion, etc. It is therefore prudent to base mass estimates of 
young stars from more robust measurements such as their spectral type, which will be done 
in a future paper. We defer further investigations of the stellar mass dependence of the disk 
fraction until this paper. 



5.3. Spatial Variations in the Disk Fraction 

The different disk fractions found in the two sub-clusters motivates an examination of 
the spatial dependence of the disk fraction. We use two methods for visualizing the spatial 
variation of the disk fraction in Cep 0B3b. First, we create maps of the disk fraction 
throughout the cluster. We created two disk fraction maps utilizing the X-ray sample and 
the visible sample. To make these maps we first created a grid covering the entire region 
with nodes spaced every ~ 50". The disk fraction at each node was then determined using 
all classified sources within a 2' radius of that node. A map was then created using the disk 
fraction at a given node as the value for the respective "pixel" on the map. Fig. [7] shows 
the results of this process for both samples, and Fig. [8] shows a 'pixel by pixel' comparison 
of the two maps. Secondly, we consider the disk fraction in annular regions in the eastern 
sub-cluster centered on the 07V star HD271086 or the western sub-cluster centered on the 
Bl V star V454 Cep, as seen in Fig. [9] for the X-ray selected set and Fig. [10] for the visible 
selected set. Note that since the background contamination may be spatially variable due to 
spatially varying extinction, we do not correct for contamination. Correcting the Class III 
sample for contamination would systematically increase the disk fraction and the amplitude 
of the observed variations. 

Both methods of visualizing the spatial variation in the disk fraction show the same 
general morphology. The eastern sub-cluster is characterized by a broad area of fairly ho- 
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mogenous disk fraction, and as shown in a recent paper by iGetman et al.l ( 120091 ). the disk 
fraction increases from about 30% to over 50% as one traces a path from the O star, to 
the Cep B molecular clump and its embedded population (Fig. [7]). The western sub-cluster, 
however, is characterized by a high disk fraction in the center of the sub-cluster that drops off 
radially outward. Note that the visible sample is biased toward the less embedded sources, 
and the map derived with this sample does not find some of the smaller clumps of stars with 
high disk fractions apparent in the X-ray map. 

One question we must ask is whether the observed variation in the disk fraction is 
consistent with random fluctuations. To address this question, we focus on the different 
disk fractions observed toward the two sub-clusters. To consider the range of disk fraction 
ratios we would expect due to statistical fluctuations, we perform a Monte-Carlo simulation. 
We assume that each member of the cluster has a probability of harboring a disk to be 
equal to the disk fraction of the entire region. For each source in a given sub-cluster we 
draw a random number from a uniform distribution between zero and one using the IDL 
routine RANDOMU. If the random value associated with a given source has a value less 
than or equal to the disk fraction for the entire sample it is considered to be a disk source, 
otherwise it is considered to be a non-disk source. Once this is done for each source, the disk 
fractions for the sub-clusters are calculated, and their ratio is determined for this particular 
realization. This is then computed for 10,000 realizations. Fig. [TT] shows the distribution of 
ratios derived using the visible sample, in the right pane, and the X-ray sample, in the left 
pane. The actual ratio is given by the red line. This exercise shows that the disk fractions 
in the two sub-cluster are different with a 95% confidence. 



5.4. Estimate of Total Cluster Membership 

We can use the disk fractions measured in the sub-clusters to estimate the total number 
of cluster members from the number of IR-excess sources. The Spitzer data shows 1135 
infrared excess sources, the vast majority of which are young stars with disks. From Fig. Q 



we asc ertain that most of the sources above 0.2 M© are detected (adopting the lBaraffe et al. 



(119981 ) tracks for a 3 Myr stars). Using the average disk fraction determined from the X-ray 
sample, 33% ± 2%, we then estimate that there are about 3000 total members in the Cep 
0B3b cluster down to the approximate mass limit of our survey. 
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DISCUSSION 



6.1. Comparison of Cep OB3b to Other Young Clusters 

Our estimated membership of Cep 0B3b of 3000 primaril y low mass stars, shows that 
Cep 0B3b is comparable in size and membership to the ONC (jAUen et al.l 120071 : ICarpenter 



2000 



1000 stars in a lOpc x 5pc region), although t he peak surface densities are con - 
siderably less dense than the peak densities in the ONC (JHillenbrand fc Hartmann I (119981 ): 
Megeath et al. (in preparation)). The content of massive stars in Cep 0B3b is only slightly 
lower tha n that of the ONC: on e 07 star in Cep 0B3b compared to an 06 and an 09.5 in 
the ONC: lMuench et al. I (120081 ). These results suggest that Cep 0B3b is the second nearest, 
large, (> 1000 star) young (< 5 Myr) stellar cluster to the Sun. Furthermore, the published 
isochrones and the relatively low extinction toward the region suggest that Cep 0B3b is 
older and more evolved than the ONC. The older ages and lower stellar ages suggest that 
this cluster has expanded since it dispersed its parental gas. 

Other than the ONC, there are no known clusters within 1 kpc that are comparable 
to Cep 0B3b in age, size and average density. The largest known young cluster next to 
the ONC and CepOB3b within 1 kpc is in the NGC 2264 region at 760 pc, part of the 
Mon OBI molecular cloud. This cluster contains ~ 450 young stars with disks, less than 



half as many as Cep03b (ISung et al.ll2009l ). The next largest known young clusters withi n 



1 kpc, ranked by number of known disk sources, are: IC 348. ~400 (jMuench et al. 



R2, 235 (IGutermuth et al.ll2009l) : IRAS 2 0050, 177 



Sigma O ri. 336 (IHernandez et al.l 120071): TR 37, ~2 00 (ISicilia-Aguilar et al 1 120061): Mon 



161 JCutermuth et al.ll2009h : IC5146, 149 JCutermuth et al.ll2009h . A major uncertainty in 



Gutermuthet al.ll2009h : AFGL 490 



2003); 



comparing the sizes of clusters is the number of diskless pre-main sequence stars. Since the 
membership counts listed in this section are primarily from Spitzer counts of disk sources, 
they are lower limits to the actual membership counts. Nevertheless, the number of known 
disk sources in these other clusters are much less than the number of known disk sources in 
Cep 0B3b. 



6.2. Disk Ablation from UV Radiation Field 



It has been known for over a decade that UV radiation from 0-stars can ablate the 
outer regions of d isks, a process that has been directly observed in the Orion Nebula Cluster 
(JBally et al.lll998l ). The IR-excesses detected by Spitzer orig inate from the warrn inner disk, 
while photoevaporation primarily truncates the outer disk (jAdams et al. II2004I ). However, 
disks may be rapidly destroyed by a combination of photoevaporation and viscous accretion 
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( iMatsuyama et al.ll2003[ ). Evidence for the destruction of inner and outer disks has been 
found in previous Spitze r observations of yo ung clusters. In NGC 2244, the Rosette Nebula, 
at a distance of 1 . 6 kpc , iBalog et al.l (|2006[ ) find outer disks being ablated near a massive O 



star. iBalog et al.l ( 120071 ) also find evidence that the fraction of young stars with disks drops 
near the rnassive O stars, evidence for the destruction of disks by the O stars. Similarly, 
Sung et al.l (120091 ) examine the young cluster NGC 2264 and find that the disk fraction tends 
to increase radially away from the 07V star S Mon. 

The only known O star within the Cep 0B3b cluster is the 07Vn star HD 217086 
(JHerrero. et al.lll992[ ). which is the primary ionizing source for the HII region Sh 155, and is 
illuminating and photoevaporating the surface of the nearby cloud as delineated by bright 
nebulosity in the IRAC and MIPS bands. HD 217086 is also situated close to the highest 
projected surface density of YSOs in the eastern sub-cluster (Fig. |3)d). A BIV star is also 
coincident with this sub-cluster. In contrast, a BOV and three BIV stars are found in the 
western sub-cluster. The Lym an continuum pro duction rate of the 07 star is nine times 
higher than that of a BOV star (jVacca et al.lll996l ). thus we must ask if the difference in the 
observed disk fractions i n the two sub-cluste rs is due to a much higher incident EUV flux in 



the eastern sub-cluster (JAdams et al. 1120041 ). 



A signature of disk photoevaporation would be a decrease in the disk fraction radially 
toward the O star. Figure [9] shows the disk fraction as a function of radial distance away 
from a massive star in each sub-cluster, using the X-ray sample. The top pane is centered on 
the 07V star HD 217086, in the eastern sub-cluster. Each annulus centered on this star has 
a width of 2' within which the disk fraction is estimated by counting the number of sources 
with and without a Spitzer excess. The bottom pane is centered on the BIV star HD 216711 
(V 454 Cep), in the western sub-cluster, and the disk fraction is estimated in the same 
manner as the top pane. The disk fraction is fairly constant as a function of radial distance 
from the 07V star in the eastern sub-cluster, while it increases significantly toward the BIV 
star in the western sub-cluster. Figure [10] shows the same radial disk fraction distribution if 
the visible sample is used to select sources. We find no evidence for disk photoevaporation 
in these plots. 

To consider whether mixing could influence the radial disk fraction distribution, i.e. 
if stars move into and away from the harsh radiation environment near the 07V star on 
timescales short enough to erase any disk fraction differences due to photoevaporation, we 
use the virial theorem to estimate the velocity of a typical star in the innermost 2' radius 
region of each sub-cluster. We assume that each star has the same mass (0.5 M©), an 
isotropic velocity distribution, and experiences the same gravitational potential. We assume 
the cluster is fully exposed, so that the mass of each sub-cluster is just the number of sources 
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times the assumed average mass and is estimated to be 46 x 0.5 M0 = 23M0 for the western 
sub-cluster and 75 x 0.5 Mq = 37.5Mq for the eastern sub-cluster. Note, here we consider 
all sources and do not apply the magnitude and color cuts of section O The typical velocity 
then would be 0.27 km s~^ in the western sub-cluster and 0.22 km s~^ in the eastern sub- 
cluster. These correspond to a time of roughly 2 Myr to move 2' in projection through the 
sub-clusters. This time is less than the lower end of the estimated age range of 3-5 Myr. 
Therefore, even though the radial disk fraction distribution near the 07V star does not 
indicate that photoevaporatation is the cause of the disk fraction difference between the two 
sub-clusters, the photoevaporation of disks in the eastern sub-cluster cannot be ruled out 
based on the radial disk fraction distribution alone. 

Finally, if disk photoevaporation is responsible for the different disk fractions in the two 
sub-clusters, we would expect the l owest mass stars to be affected the most due to their 
lower gravitational binding energy (jAdams et al. II2004I ). Instead, in Table [6l we find the 
stars with the bluest V-I colors, and hence highest masses, have the largest discrepancy in 
their disk fractions. This is evidence against disk photoevaporation. 



6.3. The Ages of Sub-Clusters and Age Spreads 



Variations in the obser ved disk fraction niay a lso be due to variations in the ages of 
the members. For instance, [Hernandez et al.l (120081 ) show that the disk fraction decreases 
from 80% to 30% as a stellar cluster/association's age increases from 1 to 5 Myr (see Fig. El) 



An age variation due to an evolutionary gradient was proposed in iGetman et al. 



Using deep IRAC observations of the eastern sub-cluster to supplement their earlier X- 
ray observations, they have evaluated disk fractions as a function of radius from the hot 
molecular core Cep B , which is in the molecular clump at the edge of the eastern sub-cluster 
(JBeuther et al.ll2000l ). They find disk fractions of ~ 75% within 3' from Cep B, ~ 55% within 
an annulus 3' to 7' from Cep B, and ~ 20% within an annulus 9' to 12' from Cep B. They 
interpret this spatial gradient as temporal, with young stars dominant in the hot molecular 
core and progressively older stars with increasing distance from the core. The temporal 
gradient results from a sequence of star formation triggered by compression driven by the 
UV radiation from the 07V star. Our eastern sub-cluster includes most of the stars in these 
annuli, and we examined the disk fraction for the X-ray sample of young stars in two sections 
within the 5' region centered on the eastern sub-cluster. The first section is toward the high 
extinction region of the cloud, and the second section is toward the lower extinction of the 
S155 HII region and HD 217086. We confirm that there is a higher disk fraction toward the 
cloud than toward HD 217086, see § 15.31 and Figures H] and [3 
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To examine the possibility that an age difference could explain the difference in d isk: 
fraction between the two sub-clusters, we use the V and I data from iMayne et al.l (120071 ) of 
the Spitzer and Chandra selected young sample. We divide the sample into two sets, one 
for each of the two 5' radius sub-cluster regions defined above. Fig. [T^ shows the location 
of these sources in V vs. V — I color-magnitude space. Note, that the reddening vector in 
V vs. V — I space is approximately parallel to the YSO locus, thus any extinction difference 
between the sub-clusters cannot masquerade as an age difference. 

By binning the V — I axis into half magnitude bins between V — I colors of 2 and 4.5, 
we find the median V magnitude of all the young sources in each bin. We then examine the 
offset of the V magnitude distribution for each sub-cluster from the median V magnitude 
by subtracting the median V magnitude of a given color bin from the V magnitude for each 
source in the same color bin. Next, we plot the normalized histogram of the distribution 
of the offsets for each sub-cluster. Due to the larger luminosity of younger stars of a given 
mass, we would expect a younger subset of stars to have a negative offset from the overall 
median. Conversely, an older subset would have a positive offset. An examination of Fig. [T2] 
shows that the sub-sample of sources from the western sub-cluster has a median offset that 
is slightly negative (-0.2 mag), while the eastern sub-clust er has an offset con sistent with 
zero. If we consider the 4 Myr isochrone from the model of iBaraffe et al.l (119981 ) and shift it 
by an amount 0.2 mag along the V magnitude axis in ^ vs. V — I color-magnitude space, 
this shift corresponds to an age difference of 1 Myr for stars within the 0.2 M© to 1.0 M© 
mass range. 



The full width half maximum of the distributions of offsets in Fig. [12] is about 1 mag 
for both sub-clusters. If we interpret this as a proxy for an age spread, then the age spreads 
for the sub-clusters are roughly 6 Myr, using the same assumptions as above. An age spread 
in the cluster is expected: we find both regions of the cluster that have dispersed their gas 
and have low disk fractions as well as regions of ongoing star formation. However, there 
are various sources of scatter in the luminosities of YSOs that could be mistaken for a n age 
spread. Potential so urces of scatter are discussed in recent papers by IJeffriej ( 120111 ) and 
Reggiani et al.l (120111 ). If this scatter is accounted for, we would expect the age spread to be 
much smaller. We will examine the observed scatter in a future paper. 



An examination of the disk fraction as a function of age curve from iHernandez et al. 



(120081 ) (see Fig. |6]) shows a roughly 13% Myr~^ drop in the disk fraction between the ages of 
2 and 6 Myr. If we assume a 1 Myr age difference between the sub-clusters, then we find a 
18% Myr~^ drop in the disk fraction using the X-ray sample and a 10% Myr~^ drop in the 
disk fraction using the visible sample. This is roughly consistent with the disk fraction curve 
from IHernandez et al.l (120081 ) . Interestingly, if we assume that there is an age difference of 
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1 Myr, then the large disk fraction difference seen in the blue bin in Table [6] could imply 
that disk dissipatio n is more rapid in the higher mass stars . This would be consistent with 



the observations of iLada et al.l ( 120061 ) and [Hernandez et al.l (120091 ) . 



7. SUMMARY 

1) We have surveyed the full spatial extent of the young Cep 0B3b cluster with the 
IRAC and MIPS instruments onboard Spitzer, and revealed over 1000 young stellar objects 
with infrared excesses. In addition we have performed Chandra imaging of the cluster with 
the ACIS instrument. Using these and archival Chandra data, we find 499 X-ray bright, 
young stars without infrared excesses. We find most of the cluster lies in a cavity of low 
extinction. The spatial distribution of these young sources shows a hierarchical morphology 
composed of two large sub-clusters with a diffuse halo. The two sub-clusters are associated 
with separate molecular clumps which contain ongoing star formation. 

2) Using these data we constrain the location of potential Cep 0B3b members in V 
and V — I color-magnitude space. We use these constraints, in addition to the known X- 
ray detections, to identify the diskless Class III sources, and conduct a census of the young 
members of Cep 0B3b. The average disk fraction for the region is 33% ± 2%. These data, 
show Cep 0B3b as a massive, nearby young cluster containing about 3000 stars, depending 
on the number of diskless pre-MS stars. 

3) Using this census of near-IR excess and diskless sources, we create a disk fraction map 
of the cluster. We find that the observed disk fraction varies spatially throughout the cluster, 
with a significant difference in the disk fractions of the two sub-clusters; 32% ± 4% in the 
eastern sub-cluster and 50% ± 6% in the western sub-cluster. This difference is dominated 
by the sources with the bluest V — I colors, corresponding to the higher mass members of 
the sub-clusters (> 1.7 M©). 

4) We compare the spatial variation in the disk fraction against the location of the most 
massive stars in the cluster and find no clear signature of disk photoevaporation. Further, 
using visible photometry, we estimate the age difference between the two sub-clusters to 
be as much as 1 Myr. The different disk fractions in the sub-clusters over this age range 
is consistent with the rate of observed disk dissipation in other star forming regions. We 
conclude that the most likely reason for the spatially varying disk fraction is a mixture of 
ages. 

5) We find that the size and membership of Cep 0B3b is similar to that of the ONC, 
making it one of the largest, known, dense, young clusters within 1 kpc of the Sun. Most of 
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the Cep 0B3b cluster lies in a low-extinction cavity, suggesting that Cep 0B3b has largely 
dispersed its natal gas and is more evolved than the more embedded ONC The size and 
evolutionary state of Cep 0B3b makes it a unique target for studying the evolution of young 
clusters. 

This work is based on observations made with the Spitzer Space Telescope, which is 
operated by the Jet Propulsion Laboratory, California Institute of Technology, under NASA 
contract 1407. Support for this work was provided by NASA through an award issued by 
JPL/Caltech. Support for this work was provided by the National Aeronautics and Space 
Administration through Chandra Award Number GO9-0017A issued by the Chandra X- 
ray Observatory Center, which is operated by the Smithsonian Astrophysical Observatory 
for and on behalf of the National Aeronautics Space Administration under contract NAS8- 
03060. STM and EK acknowledge support by NASA award 1281302. JLP acknowledges 
support from the NASA award 1276271. SJW and EDS were supported by NASA contract 
NAS8-03060 to the Chandra Science Center. Additional support for this work was provided 
by the National Science Foundation award AST-1009564. This research has made use of 
the NASA/IPAC Infrared Science Archive, which is operated by the Jet Propulsion Lab- 
oratory, California Institute of Technology, under contract with the National Aeronautics 
and Space Administration. This publication makes use of data products from the Two Mi- 
cron All Sky Survey, which is a joint project of the University of Massachusetts and the 
Infrared Processing and Analysis Center/ California Institute of Technology, funded by the 
National Aeronautics and Space Administration and the National Science Foundation and 
JPL support from SAO/JPL SV4-74011. 
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Table 1. Young Stellar Objects: Class I 



RA2000 '^ DECaooo^ V ^ V -I ^ J " H " Ks " [3.6] ^ [4.5] d [5 gj d [ggj d [34] 

22:50:14.22 62:25:50.2 ■•• ••■ •■• ■•• 15.17±0.16 14.07±0.01 13.37±0.01 12.66±0.02 11.59±0.03 8.82±0.08 0.00 

22:51:22.84 62:48:40.8 ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ 15.09±0.02 13.88±0.01 12.72±0.02 11.40±0.02 7.71±0.03 0.00 

22:52:22.88 62:28:46.6 ■■■ ■■• ■■• ■■■ •■■ 14.40±0.02 13.33±0.01 12.48±0.05 11.30±0.10 6.30±0.11 0.00 

22:52:49.10 62:29:20.5 20.89±0.01 1.96±0.02 ■■■ 15.64±0.16 14.36±0.08 12.17±0.01 11.29±0.01 10.31±0.01 9.22±0.01 6.38±0.03 0.00 

22:53:02.63 62:33:39.4 24.35±0.14 6.68±0.16 15.29±0.05 14.20±0.07 13.79±0.05 13.77±0.05 13.03±0.03 12.20±0.04 10.82±0.02 8.39±0.05 0.06 

Note. — Only a portion of the tabic is shown here. The complete table is available in the online journal. 

''RA and Dec values from the Spitzer [3.6] astrometry. 

^V and / band data from lMavne et all ||2007| '). 

"J, H and Ks band data from 2MASS. 

'^Spitzer IRAC photometry. 

'^Spitzer MIPS photometry. ^ 

'j4xa determined along the line of sight to the source using the method described in phase 2 of § [3] These values can be converted to Ay using the 
relation Ay = 9.9Axg . If ^k^ could not be estimated a value of 0.00 is shown. 

s Sources detected with Chandra are identified in the notes column with "X-ray". 
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Table 2. Young Stellar Objects: Class II 



RA2000 ^ DEC2 



yb y. 



J 



H 



Ks 



[3.6] 



[4.5] 



[5.g 



[8.0] 



24 



22 
22 
22 
22 
22 



33.73 62:48 


47.9 


36.27 62;48 


08.7 


40.12 62:56 


56.4 


47.77 62:26 


36.8 


48.63 62:26 


36.6 



15.09±0.04 14.46±0.06 14.27±0.06 14.04±0.01 13.74±0.01 14.15±0.13 • ■ ■ ■ ■ • 0.00 

12.00±0.02 11.37±0.02 10.78±0.02 9.60±0.01 8.83±0.01 8.27±0.01 7.13±0.01 4.35±0.01 0.00 

14.17±0.02 13.54±0.02 13.24±0.02 13.08±0.02 12.69±0.02 12.69±0.20 • ■ ■ ■ ■ ■ 0.00 

14.92±0.05 14.14±0.06 13.66±0.05 13.19±0.01 13.02±0.01 12.58±0.05 12.01±0.04 ■■• 0.00 

14.32±0.03 13.23±0.04 12.63±0.02 12.13±0.01 11.41±0.01 11.03±0.01 10.31±0.01 ■■• 0.04 



Note. — Only a portion of the table is shown here. The complete table is available in the online journal. 

''RA and Dec values from the Spitzer [3.6] astrometry. 

^V and I band data from lMavne et al.l |20ol). 

"J, H and Ks band data from 2MASS. 

'^Spitzer IRAC photometry. 

''Spitzer MIPS photometry. 

Aks determined along the line of sight to the source using the method described in phase 2 of §[3] These values can be converted to Ay using 
the relation Ay = 9.9j4x . If Ax could not be estimated a value of 0.00 is shown. 



to 

CO 



= Sources detected with Chandra are identified in the notes column with "X-ray". 
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Table 3. Young Stellar Objects: Transition Disk 



RA2000 " DEC 



2000 



y b V - r 



J 



H 



Ks 



[3.6] 



[4.5] 



[5.g 



[8.0] 



24 



22:48:50.00 62:33:25.7 
22:49:17.26 62:37:22.2 
22:50:44.67 62:45:28.2 
22:50:49.59 62:57:14.8 
22:50:49.95 62:31:41.1 



14.85±0.03 13.98±0.04 13.60±0.03 13.15±0.01 12.95±0.01 12.62±0.02 12.42±0.07 8.21±0.06 0.02 

10.59±0.02 10.42±0.02 10.35±0.02 10.29±0.01 10.28±0.01 10.19±0.01 9.81±0.01 7.32±0.05 0.00 

14.83±0.03 13.80±0.04 13.43±0.03 13.13±0.01 13.03±0.01 12.95±0.02 12.80±0.04 9.58±0.18 0.05 

14.64±0.04 13.76±0.04 13.48±0.03 13.18±0.01 13.01±0.01 12.83±0.03 11.78±0.04 7.36±0.03 0.03 

12.39±0.02 11.96±0.02 11.80±0.02 11.73±0.01 11.68±0.01 11.67±0.01 11.61±0.04 8.37±0.07 0.00 



Note. — Only a portion of the table is shown here. The complete table is available in the online journal. 

''RA and Dec values from the Spitzer [3.6] astrometry. 

^V and I band data from lMavne et all |20ol). 

"J, H and Ks band data from 2MASS. 

'^Spitzer IRAC photometry. 

''Spitzer MIPS photometry. 

Aks determined along the line of sight to the source using the method described in phase 2 of §[3] These values can be converted to Ay using 
the relation Ay = 9.9j4x . If Ax could not be estimated a value of 0.00 is shown. 



00 



= Sources detected with Chandra are identified in the notes column with "X-ray". 
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Table 4. Young Stellar Objects: Class III 



RA2000 '^ DECaooo^ V ^ V -I ^ J " H " Ks " [3.6] ^ [4.5] d [5 gj d [ggj d [34] 

22:51:55.39 62:41:44.7 ■•• ••■ 11.64±0.02 10.79±0.02 10.51±0.01 10.34±0.01 10.32±0.01 10.23±0.01 10.21±0.01 9.97±0.10 0.03 X-ray 

22:52:16.16 62:33:01.7 ■•• ■•■ 12.84±0.02 12.10±0.02 11.80±0.02 11.60±0.01 11.55±0.01 11.51±0.01 11.50±0.03 ■■• 0.01 X-ray 

22:52:21.1162:44:47.2 18.28±0.01 2.64±0.01 13.74±0.02 12.79±0.02 12.51±0.02 12.30±0.01 12.33±0.01 12.22±0.01 12.17±0.03 •■• 0.04 

22:52:22.12 62:34:54.3 19.99±0.01 3.26±0.01 14.31±0.04 13.22±0.04 •■■ 12.53±0.01 12.50±0.01 12.37±0.02 12.38±0.05 •■• 0.00 

22:52:22.27 62:30:33.0 18.88±0.01 2.81±0.01 13.96±0.02 12.96±0.02 12.62±0.02 12.34±0.01 12.37±0.01 12.28±0.01 12.31±0.05 ■■• 0.05 

Note. — Only a portion of the tabic is shown here. The complete table is available in the online journal. 

''RA and Dec values from the Spitzer [3.6] astrometry. 

^V and / band data from lMavne et al.l ||2007| '). 

"J, H and Ks band data from 2MASS. 

'^Spitzer IRAC photometry. 

'^Spitzer MIPS photometry. ^ 

'j4xa determined along the line of sight to the source using the method described in phase 2 of § [3] These values can be converted to Ay using the 
relation Ay = 9.9Axg . If ^k^ could not be estimated a value of 0.00 is shown. 

s Sources detected with Chandra are identified in the notes column with "X-ray". 
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Table 5. Disk Fractions in Cep 0B3b 



Method'' Total ^ Overlap " Eastern '^ Western <= Notes' 

Visible 25% ±1% 30% ±2% 35% ± 3% 45% ± 5% 

425/(425 + 1270)8389/(389 + 914)^121/(121 + 222) s«7/(87 + 106)s 
X-ray 33% ± 2% 33% ± 2% 32% ± 4% 50% ± 6% 

254/(254 + 499)g 241/(241 + 472)g 74/(74 + 151) g 80/(80 + 79)E 

Visible 32% ± 2% 36% ± 2% 39% ± 4% 53% ± 6% Sub 

425/(425 + 889)g 389/(389 + 914)gl21/(121 + 191) g87/(87 + 75)g 
X-ray 36% ± 2% 36% ± 2% 34% ± 4% 52% ± 6% Sub 

254/(254 + 446)g 241/(241 -|- 431)g 74/(74 + 145) g 80/(80 + 73)g 

Visible 24% ± 1% 29% ± 2% 34% ± 3% 44% ± 5% Cut 

367/(367 + 1139)g338/(338 -I- 811)gl03/(103 -I- 199)g 75/(75 -I- 96)g 
X-ray 33% ± 2% 33% ± 2% 32% ± 4%^ 50% ± 6%^ Cut 

223/(223 + 451)g 214/(214 -|- 426)g 64/(64 + 137)g 71/(71 + 71) g 

''Method used to determine potential Class III members. 

''The full field of view for the given sample. 

"^The field of view common to both the visible and the X-ray sample. 

'^5' radius regions centered on RA; 22:56:49.44, Dec: 62:39:55.54. 

°5' radius regions centered on RA: 22:53:50.83, Dec: 62:35:36.37. 

Samples with estimated contaminants subtracted are identified with the note Sub. 
See § 13.21 for details about the X-ray contamination and § 13.31 for details about the 
visible contamination. Samples with the NIR color and magnitude cut applied are 
identified with the note Cut. See § 15.11 for details of the NIR color and magnitude 
cuts. 

gThese numbers are (number of disk sources) /(number of disk sources + number of 
non-disk sources) 



''These are the disk fractions quoted in the abstract. 
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Table 6. Disk Fractions in the Sub-Clusters as a Function of Dereddened V — I color 

Color Range ^ Spectral Type'' Mass '^ Sample '^ Overlap ° Eastern Western ^ 

-0.5 to 1.05 mag > K3 > 1.7 Visible 19% ± 2% 19% ± 4% 39% ± 9% 

72/(72 + 311)'^ 18/(18 + 78)''20/(20 + 31)^ 
Visible Subtracted 30% ± 4% 23% ± 6% 63% ± 14% 

72/(72 + 172)'' 18/(18 + 59)''20/(20 + 12)* ^ 

X-ray 30% ± 5% 15% ± 6% 56% ± 14% 

44/(44 -f 103) ^ 8/(8 + 43)'' 15/(15 + 12)' ' 

1.05 to 1.72 mag K3 to MO 1.7 to 0.6 Visible 31% ± 3% 40% ± 6% 39% ± 9% 

120/(120 -I- 263)'i40/(40 + 59)'25/(25 -I- 39)'' 
Visible Subtracted 38% ± 3% 44% ± 7% 45% ± 9% 

120/(120 -I- 199)''40/(40 + 50)''25/(25 -I- 30)' ' 

X^ 32% ± 4% 33% ± 7% 46% ±11% 

58/(58 -I- 124)'' 21/(21 + 42)''l 7/(17 -I- 22)'' 
1.72 to 4.35 mag MO to M6 0.6 to 0.2 Visible 38% ± 3% 42% ± 6% 54% ± 10% 

146/(146 -I- 237)''45/(45 + 62)''30/(30 + 26)'' 
Visible Subtracted 41% ± 3% 43% ± 7% 57% ± 10% 

146/(146 -I- 237)''45/(45 + 59)''30/(30 + 23)'' 
X-ray 31% ± 5% 26% ± 8% 50% ± 16% 

36/(36 -I- 80)'' 10/(10 -I- 28)''10/(10 -I- 10)'' 



Note. — Due to the unknown V — / color dependence of the estimated X-ray contamination, only the 
visible sample has contamination subtracted. 

"The total sample was divided into blue, intermediate and red thirds in the dereddened V — I color space. 
The colors given are the bounding colors. Also given is the corresponding spectral type and corresponding 
mass range for these bounding colors. See ij 15.21 for details. 

''Spectral types correspond to the bounding colors are from iKenvon Sz Hartmann I 1 119951) . 

"Masses corresponding to the bounding colors are estimated fromthe models of lSiess et al.l 1 120001 ). 

''Sample used to determine potential Class III members. 

■^AU sources in the field of view common to both samples. 

''5' radius region centered on: RA: 22:56:49.44, Dec: 62:39:55.54. 

"5' radius region centered on: RA: 22:53:50.83, Dec: 62:35:36.37. 

''These numbers are (number of disk sources) /(number of disk sources -|- number of non-disk sources) 
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22" 58"'5 



Fig. 1. — Top: Three-color 3.6 (blue), 5.8 (green), and 24 /im (red) image of Ceplieus 0B3b 
region. The extended pink emission is likely a mixture of emission from PAHs (particularly 
the 6.2 fim feature) mixed with thermal emission from dust grains in the 24 /xm band. This 
emission defines the border of a cavity in the molecular cloud within which the young cluster 
is located. The extended red emission is due to warm dust at 24 /im in the HIT region. The 
massive stars in the region are labeled; the 07V star is HD 217086 and the early B stars 
are HD 217061, BHJ 22, BHJ 24, BHJ 17br and V454 Cep. Bottom: Color composite Poss 
II image showing Cep 0B3b in the visible. The nebulous green emis sion is the HII region 

( 2006) 
fl2007h 



Sharpless 155. The regions imaged by Chandra for this study and from lCetman et al. 



are outlined with the solid white line. The region imaged in the visible by lMayne et al. 
is outlined with the dashed white line. 
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Fig. 2. — IRAC/MIPS color-color and color-magnitude diagrams of Cep 0B3b. The colors of 
the markers correspond to the classifications determined with the three phases from Sec. 3. 
Red denotes Class I objects, green denotes Class II objects, blue denotes transition disk 
objects and black denotes Class Ill/field stars. The dashed line in the [4.5]- [8] vs. [4.5] 
diagram (lower left plot) separates the region occupied by AGN contaminants. All objects 
in this region were flagged as contaminants and deleted. They were later included only if they 
have bright MIPS 24/xm detection ([24] < 7) or very red IRAC/MIPS colors ([3.6] -[5.8] > 0.5 
and [4.5] — [24] > 4.5 and [8.0] — [24] > 4). The black arrows show the reddening vectors for 
an Ak = 5. 
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Fig. 3. — V Ys. V — I color-magnitude diagram of Cep 0B3b. The black dots are sources in 
the Cep 0B3b field with photometry that meets the photometric error threshold described in 
the text, and the blue dots are sources classified as pre-MS. The region bounded by the black 
curve represents the YSO isochrone as described in the text. The dashed white and black 
line are the ZAMS of ISiess et al.l (120001). and the black a r row c orresponds to 1 magnitude of 
extinction using the reddening law of iRieke fc LebofskyI (Il985l ). The left plot shows sources 
classified based using Spitzer colors and the right plot shows sources classified based on their 
detection in X-rays. 
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Fig. 4. — Top: grayscale extinction map of same region with countours corresponding to 
Ay of 2.5, 5.0, 7.5, and 10. The red and green dots are Spitzer classified young stars with 
circumstellar excess indicative of an envelope and a disk, respectively, and the blue dots are 
Chandra detected young stars with increased coronal activity. The large red circles show the 
adopted boundaries of the sub-clusters. Bottom: YSO surface density map created using 
IR, X-ray and visible selected YSO candidates. YSO surface density was calculated using 
the 11 nearest neighbors. Contours correspond to YSO surface densities of 10, 30 and 70 
stars pc~^. 
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Fig. 5. — J vs. J — H color-magnitude diagram of Cep 0B3b. The black dots are sources in 
the Cep 0B3b field with photometry that meets the photometric error threshold described 
in the text, while the green, red and blue dots are further selected as young sources based 
on Spitzer, Chandra and/or visible photometry. The red and green sources are excess and 
non-excess sources, respectively, that satisfy the magnitude-extinction requirements for the 
disk fraction analysis in section |5l while the blue sources do not. The upward sloping dashed 
lines represent Ay values of 0, 2, 4, 6, and 8 magnitudes. The downward sloping solid lines 
are reddening vectors for a mj = 13.01, 13.94, 14.75, and 15.81 star. These magnitudes 
correspond to a 0.7, 0.35, 0.2 and 0.08 Mp-, s t ar at an age of 3 Myr at the adopted distance 
of 700 pc using the model of iBaraffe et al.l (Il998[ ). Each pane shows a different selection 
method. From left to right: Spitzer infrared excess sources, Chandra X-ray sources and 
visible wavelength selected sources. The criteria for each selection method are described in 
the text. 
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Fig. 6. — Disk fractio n vs. age for young stell ar clusters in the literature. This plot is taken 
from Figure 14 from iHernandez et al.l (120081 ) and augmented with the disk fractions and 
assumed age of the Cep 0B3b sub-clusters in red and the average disk fraction and assumed 
age of the entire Cep 0B3b region in green. 
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Fig. 7. — Top: Disk fraction map created with the X-ray selected members. Contours 
represent disk fractions of 20, 34, and 50 percent. Red circles outline the sub-clusters and 
the white dots are the locations of the massive stars. Bottom: Disk fraction map created 
with the visible selected members. Contours are at the same levels as above. 
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Fig. 8. — A pixel by pixel comparison of the disk fraction maps created with Chandra selected 
YSOs and visible selected YSOs. The pixels considered are those with a measurement 
uncertainty of less than 10 percent. The dashed line, with a slope of unity, is plotted for 
reference. 
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Fig. 9. — Disk fraction as a function of radial position for the X-ray selected young stars with 
(filled red circles) and without (open blue circles) disks, subject to the mass and extinction 
requirements discussed in the text. The top plots are centered on the 07V star HD 217086, 
while the bottom lots are centered on V454 Cep, a BIV star at the center of the western 
sub-cluster. The disk fraction is determined within each annuli, which are separated by 2'. 
Note, the disk fraction in the outer annulus of the western sub-cluster is increased by the 
small grouping of stars with disks between the two sub-clusters. 
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Fig. 10. — Disk fraction as a function of radial position for the visible selected young stars 
with (filled red circles) and without (open blue circles) disks, subject to the mass and ex- 
tinction requirements discussed in the text. The center positions and annuli are the same as 
Figure O Note, that the visible sample does not pick up the same increase in disk fraction in 
the outer annuli as the X-ray sample. This could be due to preferentially selecting against 
reddened disk sources when requiring a V^-band detection. 
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Fig. 11. — Results of a Monte Carlo simulation of the disk fractions in the sub-clusters. The 
histogram show the results for 10,000 realizations of the simulation, the solid black line is 
the mean, the dashed black lines are 2cr from the mean and the red line is the actual ratio. 
The left pane shows the X-ray selected sources, and the right pane shows the visible selected 
sources. 
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Fig. 12. — Left: l^ vs. V — I color-magnitude diagram for the members of the eastern 
sub-cluster, open circles, and western sub-cluster, filled circles. Filled red circles are the 
median V magnitudes for the color bins and the red lines show the width of each bin. 
Right: distribution of V magnitude deviations from the median for the young sources in 0.5 
magnitude color bins. The dashed and solid red lines show the median color deviation for 
the eastern and western sub-clusters, respectively. 



